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Numerical  analysis  of  near-field  optical  trapping 
using  tapered  fiber  probe 

Guoping  Zhang*,  Zhongru  Zhu,  Yanping  Li,  Ge  Xia,  and  Qi  Lin 
Department  of  Physics,  Central  China  Normal  University,  Wuhan,  430079,  China 


ABSTRACT 

A new  scheme  for  optical  trapping  is  presented  in  this  paper.  The  method  is  based  on  a tapered  fiber  probe  with  a tip 
diameter  less  than  a light  wavelength.  A three-dimension  gradient  optical  field  is  formed  within  the  optical  near  field  of 
the  fiber  probe,  and  a particle  approaching  the  fiber  probe  tip  will  be  trapped.  The  evanescent  electromagnetic  field  in 
the  vicinity  of  the  fiber  tip  is  calculated  by  the  multiple  multipole  method  (MMP).  The  intensity  distributions  and  the 
trapping  potential  of  the  near  fields  of  the  tip  versus  the  longitudinal  and  transverse  distances  from  the  tip  are  analyzed 
respectively.  The  trapping  force  is  obtained  for  a dielectric  particle.  The  numerical  calculating  results  show  the 
availability  of  this  method. 

Key  words:  optical  trap,  near-field  optics,  fiber  probe,  scanning  near-field  optical  microscope 

1.  INTRODUCTION 

Optical  trapping  by  a single-beam  gradient  force  trap  was  proposed  by  Ashkin  et  al.  for  the  first  time.111  Since  then,  this 
method  has  been  developed  as  an  optical  tweezers  technology,  and  extensively  used  for  manipulation  of  various 
submicron-size  dielectric  particles  and  biological  structures  (such  as  viruses,  bacteria  and  yeast  cells).12,31  However, 
conventional  optical  tweezers,  based  on  the  gradient  field  near  the  focus  of  a laser  beam,  have  a diffraction-limited 
trapping  volume. 

Scanning  near-field  optical  microscopy  (SNOM),  developed  recently  from  the  combination  of  the  scanning  probe 
technology  with  the  optical  microscopy,  offers  an  optical  image  technology  with  an  ultra-high  resolution  beyond  the 
diffraction  limit.14'61  The  resolution  capability  of  the  SNOM  system  is  determined  by  the  transmission  property  of  the 
fiber  probe,  such  as  the  output  photon  flux  density  from  the  probe  and  the  tip  size  of  the  probe.  Nowadays  its 
longitudinal  and  transversal  resolutions  have  reached  O.lnm  and  7nm  respectively. 

In  this  paper,  a novel  high-resolution  optical  trapping  scheme  is  presented  by  a tapered  fiber  probe  used  for  SNOM. 
Because  the  tip  diameter  of  the  fiber  probe  is  smaller  than  a light  wavelength,  the  near  field  close  to  the  tip  mainly 
consists  of  evanescent  components  which  decay  rapidly  with  distance  from  the  tip.[7'91  Thus  a three-dimension  gradient 
optical  field  is  formed  within  the  optical  near  field  of  the  fiber  probe.  When  one  particle  approaches  the  fiber  probe  tip,  it 
will  be  trapped  by  a three-dimension  optical  trap,  which  is  produced  from  the  action  of  the  three-dimension  optical 
gradient  force.  As  this  type  of  optical  trapping  appears  within  the  optical  near  field  of  the  fiber  probe  tip,  it  is  called  as 
the  near-field  optical  trap. 
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2.  TAPERED  FIBER  PROBE 


Tapered  fiber  probes  used  for  the  scanning  near-field  optical  microscopy  are  generally  fabricated  by  local  melting  and 
pulling  single-mode  optical  fibers.110,1'5  An  about  2-meter-length  single-mode  fiber  stripped  off  a short  section  of  jacket 
in  its  middle  part  is  heated  locally  by  a beam  of  C02  laser,  simultaneously  an  external  pulling  force  is  applied  on  two 
ends  of  the  fiber.  When  it  is  broken,  a tapered  fiber  probe  is  obtained.  Appropriate  choosing  the  laser  intensity,  the 
heating  length  and  time,  and  the  pulling  force  strength  can  be  used  to  yield  fiber  probes  with  various  tip  sizes  and  taper 
shapes.  The  tip  diameters  of  the  fiber  probes  used  for  the  SNOM  system  are  required  less  than  an  optical  wavelength, 
e.g.,  tens  to  hundreds  of  nanometers  for  the  visible  light.  Figure  1 gives  schematic  diagrams  of  the  conical-taper  and  the 
parabolic-taper  fiber  probes. 


Fig.l.  Schematic  diagram  of  tapered  fiber  probes  with  (a)  a conical  taper  and  (b)  a parabolic  taper.  The 
areas  within  the  dash  lines  stand  for  the  fiber  core,  and  those  between  the  solid  lines  and  the  dash  lines  for 
the  fiber  cladding. 

Both  the  diameters  of  the  cladding  and  the  core  of  the  tapered  fiber  probe  decrease  gradually  along  the  taper,  and  they 
are  usually  assumed  to  maintain  their  initial  ratio.  The  modes  in  the  beginning  of  the  fiber  taper  exhibit  first  the  core- 
guided  modes,  and  they  contracts  slightly.  After  propagating  a certain  distance  along  the  taper,  the  core-guided  modes 
spread  into  the  cladding,  and  gradually  couples  to  the  cladding-guided  modes,  for  which  the  boundary  condition 
becomes  the  cladding/air  boundary. 

The  waveguide  structure  of  the  tapered  fiber  is  changed  with  a different  taper  shape,  thus  results  in  different  transmission 
property.  It  is  shown  by  the  local  modal  analysis  that,  the  mode  field  is  first  the  core-guided  mode  and  contracts  slightly 
along  the  taper,  then  spreads  into  the  cladding  and  gradually  couples  to  the  cladding-guided  modes.  When  the  tip  size  of 
the  tapered  fiber  is  less  than  an  optical  wavelength,  the  evanescent  wave  becomes  the  major  part  of  the  optical  near  field 
from  the  tapered  fiber.  Meanwhile,  due  to  the  tip  scattering,  the  evanescent  wave  is  mostly  concentrated  at  the  axial 
direction,  and  the  lateral  optical  field  decreases  gradually.  So  optical  field  gradients  are  formed  near  the  tip  of  the 


322 


tapered  fiber  probe. 

The  transmission  properties  and  the  electromagnetic  field  distributions  are  different  for  the  conical  and  parabolic  tapered 
fiber  probes.191  Generally,  the  parabolic  tapered  probe  has  a transmission  efficiency  of  almost  one  order  of  magnitude 
higher  than  the  conical  one.  The  electromagnetic  near-field  distributions  of  these  two  probes  are  more  focused  while 
their  taper  angles  increase. 


3.  THEORY 

A laser  beam  is  coupled  into  the  single-mode  fiber  from  the  end  opposite  to  the  tip.  In  order  to  solve  Maxwell’s 
equations  in  the  specific  geometry  of  the  tip  of  the  tapered  fiber  probe  and  its  environment,  we  employ  the  multiple 
multipole  method  (MMP)  which  recently  has  been  applied  to  study  near-field  optical  phenomena.112,131  The  MMP 
method  expands  the  solutions  on  multipolar  eigenfunctions.  The  vector  multipolar  eigenfunctions  F„(r,®)  for  the 
eigenvalue  q„  satisfy  the  vector  wave  equation: 

-VxVxF„(r,ty)  + ^2F„(r,«)  = 0.  (1) 

In  terms  of  scalar  fields,  the  scalar  multipolar  eigenfunctions  $,(r,m)  for  the  eigenvalues  q„  satisfy  the  scalar  Helmholtz 
equation: 

VV„(r,®)  + ?„V»(r>®)  = 0,  (2) 

where 

jV,*  (r,  of)  </>„■  (r  ,a)dr  = Snn.  (3) 

These  eigenfunctions  form  an  orthonormal  basis  set  in  the  Hilbert  space.  The  simplest  solution  is  obtained  in  Cartesian 
coordinates: 

4,n(r >a)  = ~7=~  exPOk ' r}  > (4) 

V 8;r3 

where  k is  the  wave  vector.  The  multipolar  wave  functions  are  formulated  in  the  spherical  coordinates: 

w ^ fcOS/WC?] 

0n(r,a))  = <tal  (r,co)  = P,  (cos d)z,(qnr)\  . (5) 

[sin/w^J 

where  P”  (cos  0)  represents  for  the  associated  Legendre  polynomials  and  zi(q„r)  for  the  spherical  Bessel  functions.  The 
index  cr distinguishes  between  the  even  and  odd  functions.  In  the  cylindrical  coordinates,  the  multipolar  wave  functions 
are  characterized  as 

4>n(r,co)  = <j>a  m qn  (r,co)  = B„(k  p)exp{jq„zz}\^ 
where  the  Bm  are  Bessel  functions,  k = Jq%  x + ql  y and  p = yjx2  + y2  . 

Three  sets  of  vector  eigenfunctions  can  be  built  by  applying  the  gradient  operator  to  the  scalar  functions: 


cos  m<p 1 
smm(p  \ ’ 


(6) 
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L„(r,®)  = CV$i„(r,®), 

(7) 

M„(r,®)  = C V x^„(r,®)S , 

(8) 

N„(r,®)  = 7VxVx^(r,a))S, 
k 

(9) 

where  C is  a normalization  factor  which  relies  on  the  coordinates  system,  and  S is  a constant  vector  of  unit  length, 
sometimes  called  “Poynting  vector”. 

Because  of  the  property  of  the  Poynting  vector  and  the  orthonormalization  of  the  scalar  eigenfunctions,  the  above  three 
sets  L„(r ,co),  M„(r ,co)  and  N„(r,a>)  are  mutually  orthogonal  in  the  Hilbert  sense.  The  following  relationship  can  be  easily 
proved  valid  for  an  infinite  homogeneous  system: 

[l„  (r,  co ) • L*„  (r\ co)  + M„(r,»)  • M*  (r',  co)  + N„(r,  co)  • N*  (r%  ®)]  =£(r  - r') . (10) 

n 

The  longitudinal  eigenfunctions  L„(r,®)  have  no  physical  meaning  for  our  discussion.  Therefore,  after  dividing  the  space 
in  subdomains  where  the  refractive  index  is  constant,  the  MMP  method  performs  the  expansion  of  the  electric  field  in 
each  subdomain  only  on  the  sets  of  transverse  eigenfunctions  M„(r,®)  and  N„(r,®): 

E'(r ,co)  = ^KPMn(r,®)  + ^N„(r,®)] . (11) 

n 

The  unknown  coefficients  an  and  b„  are  then  determined  from  the  electromagnetic  boundary  conditions  on  the  interfaces 
between  adjacent  subdomains  by  a least-squares  minimization. 

The  multipolar  functions  used  in  the  basis  set  of  the  multiple  multipole  method  are  so  short  range  that  only  their  close 
neighborhood  is  affected.  Thus  the  method  is  better  suited  to  account  for  localized  geometries  than  the  expansion  in 
plane  waves. 

While  the  electric  field  distribution  near  the  tip  is  determined,  the  gradient  force  for  a Rayleigh  particle  can  be  easily 
calculated  as 

F = |V£2,  (12) 

where  a is  the  polarizability  of  the  particle.  The  particle  tends  to  move  to  the  higher  intensity  region  where  its  induced 
dipole  has  lower  potential  energy.  To  overcome  the  limitation  of  Rayleigh  approximation,  a rigorous  treatment  of  the 
trapping  force  can  be  performed  by  applying  the  conservation  law  for  momentum. 

The  gradient  force  can  be  expressed  by  the  Maxwell’s  stress  tensor  T as 

F = ^ < T ■ n >dS , (13) 

where  <...>  stands  for  the  time  average,  ^denotes  a surface  enclosing  the  particle,  n is  the  outwardly  directed  normal 
unit  vector,  and 

T = £-0£-EE  + //oPHH  - ^(£0££2  + • (14) 
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Where,  I stands  for  the  unit  dyad  and  s,  /i  are  the  dielectric  constant  and  magnetic  permeability  of  the  surrounding 
medium,  respectively.  For  a particle  in  the  vicinity  of  the  tip,  the  magnetic  contribution  to  the  force  is  found 
approximately  two  orders  of  magnitude  lower  than  the  electric  one. 

The  trapping  potential  U of  a particle  located  at  r0  is  then  given  by 

U(r0)  = -f°F(r)-</r.  (15) 

J oo 


4.  RESULTS  AND  DISCUSSIONS 

The  following  parameters  of  the  single-mode  fiber  are  chosen  in  our  calculation.  The  refractive  indices  of  the  core  and 
cladding  are  ni=1.4544  and  n2=1.45,  respectively.  The  diameters  of  the  core  and  cladding  are  respectively  2a=8pm  and 
2fr=125pm.  The  wavelength  of  the  illuminating  light  is  1=632. 8nm  (He-Ne  laser).  The  fiber  taper  is  assumed  to  be  linear 
with  an  lOOnm-tip  size. 

Figure  2 shows  the  contours  of  the  near  field  distribution  in  the  vicinity  of  the  tip  of  the  fiber  probe  with  a conical  taper. 
The  near  field  close  to  the  tip  is  mainly  of  evanescent  components  which  are  mostly  concentrated  at  the  axial  direction. 
The  field  is  rotationally  symmetric  in  the  vicinity  of  the  tip,  and  attenuate  rapidly  with  distance  from  the  tip.  So  a three 
dimensional  gradient  optical  field  is  formed  near  the  tip  of  the  tapered  fiber  probe.  The  field  is  rotationally  symmetric  in 
the  vicinity  of  the  tip  since  the  fiber  probe  has  a rotationally  symmetric  geometry. 


Fig.2.  Contours  of  the  near  field  of  tapered  fiber  probe.  The  single-mode  fiber  is  specified  as  n,=1.4544,  n2=1.45,  2a=8pm,  and 
26=125pm.  The  illuminating  wavelength  is  1=632. 8nm  (He-Ne  laser).  When  a particle  approaches  the  vicinity  of  the  fiber  tip,  it 
will  be  trapped  by  the  gradient  force  of  the  near  field  of  the  fiber  probe. 

Figures  3 and  4 give  the  intensity  distributions  of  the  near  fields  out  of  the  tip  as  functions  of  the  longitudinal  and 
transverse  distances  (z  and  x,  respectively)  from  the  tip,  respectively.  Since  the  optical  near  field  in  the  vicinity  of  the  tip 
is  evanescent,  the  field  amplitude  reduces  dramatically  with  the  distance  from  the  tip,  and  shows  almost  rotationally 
symmetric  in  the  vicinity  of  the  tip. 
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Fig.3.  Intensity  distribution  of  the  near  field  of  the  tip  along 
the  longitudinal  direction  z.  The  illuminating  wavelength  is  A 
= 632.8nm  (He-Ne  laser).  The  tip  diameter  of  the  fiber  probe 
is  assumed  as  lOOnm. 


Fig.4.  Intensity  distribution  of  the  near  field  of  the  tip  along 
the  transverse  direction  x with  different  longitudinal  distances 
(z=0,  lOnm  and  20nm,  respectively).  The  field  is  rotationally 
symmetric  in  the  vicinity  of  the  tip. 


In  our  calculation,  the  interaction  of  the  tip  and  a nanometric  particle  close  to  the  tip  is  not  counted.  However,  the  field 
distributions  will  be  distorted  around  a dielectric  particle.  A rigorous  treatment  of  the  gradient  field  and  the  trapping  force 
may  be  performed  by  applying  the  conservation  law  for  momentum.1’41 

When  a particle  (of  diameter  rf=T0nm)  approaches  the  fiber  tip,  the  trapping  force  and  the  potential  energy  can  be 
calculated  from  the  above  equations  (12)  to  (15).  Figures  5 and  6 show  cross  sections  of  the  trapping  potential  along  the 
axial  (z)  and  lateral  (x)  directions  from  the  tip,  respectively.  Since  the  trapping  potential  is  almost  rotationally  symmetric, 
we  only  show  the  results  for  the  (x,  z)  plane. 


Fig.5.  Normalized  trapping  potential  energy  of  a particle 
(of  diameter  d=10nm)  as  the  function  of  the  longitudinal 
distance  z from  the  tip. 


Fig.6.  Normalized  trapping  potential  energy  of  a particle 
(rf=10nm)  as  the  function  of  the  transverse  distance  x from 
the  tip,  with  different  longitudinal  distances  (z=0,  lOnm  and 
20nm,  respectively). 
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The  trapping  potential  is  sensitively  dependent  upon  the  size,  shape  and  dielectric  constant  of  the  trapping  particle,  as  well 
as  the  size  and  refractive  index  of  the  fiber  tip.  Generally,  the  trapping  force  is  inversely  proportional  to  the  tip  size,  and  a 
sharper  tip  is  required  for  trapping  smaller  particles. 

The  trap  force  may  be  affected  by  the  wavelength  of  the  incident  light.  The  particle  can  be  trapped  and  approach  the  probe 
for  one  wavelength,  however,  it  is  also  possible  to  be  repelled  and  escape  from  the  probe  for  another  wavelength.  By  the 
wavelength-dependent  optical  force,  one  can  trap,  move,  and  deposit  the  particle  on  a desired  position  on  a nanometer 
scale. 


5.  CONCLUSIONS 

A method  performing  near-field  optical  trapping  by  a tapered  fiber  probe  has  been  presented  in  this  paper.  Since  the  fiber 
probe  has  a nanometer-scaled  resolution,  the  near-field  optical  trap  is  better  suited  to  trap  smaller  particles.  The  multiple 
multipole  method  is  employed  to  calculate  the  near  field  in  the  vicinity  of  the  fiber  tip  and  obtain  the  trapping  force  for  a 
dielectric  particle.  The  calculation  results  indicate  the  feasibility  of  the  method. 

Compared  with  the  optical  trap  by  a highly  focused  laser  beam,  the  near-field  optical  trap  using  a tapered  fiber  probe 
exhibits  the  following  advantages:  (1)  the  optical  trapping  system  is  simplified  without  focusing  objective  lenses;  (2)  the 
fiber  is  conveniently  coupled  with  a laser  diode  acted  as  the  illuminating  source,  so  that  high  a coupling  efficiency  is  easily 
obtained;  (3)  the  highly  confined  evanescent  fields  significantly  reduce  the  trapping  volume;  (4)  the  large  field  gradients 
result  in  a large  trapping  force;  and  (5)  the  trapped  particles  can  be  moved  precisely  and  freely  when  the  fiber  probe  is 
mounted  on  the  three  dimensional  scanning  set-up  of  the  SNOM  system. 

The  experiment  of  the  near-field  optical  trapping  is  in  progress,  and  the  experimental  results  will  be  reported  in  another 
article  in  the  near  future. 
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